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ABSTRACT 


Television and infrared sensing techniques utilized by meteorological satellites pro- 
vide a very useful tool for oceanographic research. 

Satellite television pictures of cloud patterns overlying thermally different oceanic 
boundaries reveal significant features of sea surface temperature distribution. Nimbus II 
High Resolution Infrared Radiometer (HRIR) data which is sensitive in the 3.4 to 4.2 
micron window were analyzed over several oceanic regions. Quantitative sea surface 
temperatures within an accuracy of ±2°K were derived from nighttime satellite radiance 
measurements in the absence of clouds. Comparisons with observations from low flying 
radiometer-equipped aircraft confirmed good absolute sea surface temperature measure- 
ments from the satellite radiometer if 3° to 4°K corrections for atmospheric attenuation 
and other losses were applied. Sea surface temperature discontinuities such as the North 
Wall of the Gulf Stream have been located consistently within 10 kilometers of the posi- 
tions indicated by the airplane radiation data. With present techniques, primarily de- 
signed for meteorological purposes, the Gulf Stream boundary has been seen, at least in 
significant parts in about 50 out of 175 days. Similar results have also been obtained in 
analyses of the Agulhas Current boundary, the boundary between the Brazil and Falkland 
Currents and the Kurishio - Oyashio mixing region. These satellite radiation observa- 
tions suggest that the Brazil- Falkland Current boundary which is associated with a sur- 
face temperature discontinuity is as sharp and strong as the Gulf Stream North Wall. The 
Agulhas Current exhibits a similar temperature discontinuity along its western boundary, 
which separates it from the Benguela Current surface waters. 

Daytime observations within the 3.4 to 4.2 micron window have also shown qualita- 
tively the location of major current boundaries. 

For detailed oceanographic research work in the future, a higher spatial resolution 
and a shift to the 10 to 11 micron region is desirable. Simultaneous window channel ob- 
servations with a visible channel (0.55 to 0.75 m) would be necessary to delineate and 
eliminate cloud contaminated areas within the field of view of the radiometer. 
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OCEAN CURRENT AND SEA SURFACE TEMPERATURE 
OBSERVATIONS FROM METEOROLOGICAL SATELLITES 


by 

Guenter Warnecke,* Larry M. McMillin,* and Lewis J. Allison 
Goddard Space Flight Center 


INTRODUCTION 

Polar -orbiting and geostationary satellites are valuable observational tools in geophysics be- 
cause of their ability to scan the distribution of physical properties in various scales synoptically 
over large areas or over the entire globe in a reasonably short time. In probing the earths sur- 
face and lower atmosphere from space, the remote sensing of electromagnetic radiation is of pri- 
mary importance. The limited penetration of electromagnetic waves through water restricts 
oceanographic applications to relatively shallow layers or even to the immediate surface only, de- 
pending on the spectral region. 

Although no specific oceanographic experiments have yet been performed from spacecraft, 
airborne oceanographic observations from airplanes as well as from manned spaceflights and un- 
manned meteorological satellites indicate a large potential for future oceanographic research and 
application. In the visible part of the spectrum, visual and photographic observations from air- 
planes and manned spacecraft have shown evidence that remote aerial observations can contribute 
considerably to near -coast and open-sea research (Ewing, 1965; National Council on Marine Re- 
sources and Engineering Development, 1967). A number of experiments have been performed in 
the infrared to measure, for meteorological purposes, the radiation of the earth-atmosphere sys- 
tem in various spectral regions (Bandeen, 1968). In an atmospheric "window,” i.e., in a spectral 
region where gaseous absorption does not exist (ideal case) or is at a pronounced minimum 
(natural case), according to radiative transfer theory, measurements of the outgoing thermal radi- 
ance can be interpreted in terms of the temperature of an opaque surface (e.g. land, water, thick 
clouds, etc.) provided that the emissivity of the surface is known and the surface fills the field of 
view of the radiometer. Therefore, under clear-sky conditions, radiometric measurements in a 
spectral window region taken at satellite altitudes should provide quantitative information on sea 
surface temperature. Because of the quasipolar orbits of present weather satellites and the 
scanning mode of the radiometers used, global mapping of sea surface temperatures should be- 
come possible. 


*On leave from the Freie Universitat Berlin, Germany, as a NAS-NRC Senior Post-doctoral Resident Research Associate with the 
Planetary Radiations Branch, Goddard Space Flight Center. 

tAllied Research Inc., Concord, Mass.; present affiliation, Department of Earth Sciences, Iowa State University, Ames, Iowa. 
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Infrared techniques in measuring sea surface temperature have been applied before from ship- 
board and airplanes (J. Clark, 1964; Wilkerson, 1966), and a detailed discussion of the physics 
of aircraft measurements has been given by H. L. Clark (1967). However, only sporadic measure- 
ments are possible from a ship, and even airborne scanning radiometers can only cover very 
limited areas that are restricted by the limited flight altitude, operational range of the aircraft, 
and frequency of flights. Thus, a polar- orbiting satellite with its ability of repeatedly mapping the 
entire globe in 12-hour intervals suggests an ideal means of monitoring sea surface temperature 
distributions from space. 

In this report, the remote detection of ocean currents and sea surface temperature patterns 
will be discussed. A brief description of the theoretical background of sea surface temperature 
measurements from space will be given, and the feasibility of this method will be explored by 
means of experimental results mainly from measurements aboard the Nimbus II satellite. 


INFERENCE OF SEA SURFACE TEMPERATURE PATTERNS FROM 
SATELLITE CLOUD OBSERVATIONS 

Slight color differences along the boundary between ocean currents have been observed from 
aircraft and documented in color photographs; however, no multispectral color techniques have 
been developed for spacecraft utilization. Present photographic (television) techniques of obser- 
vation from satellites designed merely for meteorological purposes do not have the spatial and 
spectral resolution necessary for such a purpose. Thus, at present only inferences of differential 
sea surface temperature from televised cloud features are possible. This method is, however, 
restricted to certain meteorological conditions. 

The cloud patterns seen in television pictures are determined mainly by the structural and 
dynamical conditions within the atmosphere itself. The influence of surface properties on the cloud 
structure is extremely variable in character, space, and time. It can be recognized easily in the 
case of temperature differences between land and the oceans. However, over the ocean, with its 
comparably small horizontal temperature differences, the surface effect usually appears over- 
ridden by the effect of internal atmospheric processes on the cloud structure. In the case of rather 
homogeneous horizontal conditions and fairly weak large-scale dynamic activity in the atmosphere, 
sea surface temperature characteristics can be reflected in the cloud patterns, as has been shown 
by Warnecke et al. (1968) and the National Council on Marine Resources and Engineering Develop- 
ment (1967). An excellent example is presented in Figure 1. 

On April 10, 1967, the satellite ATS-II made a remarkable TV photograph of the eastern tropical 
Pacific Ocean. In this photograph, the position of the Humboldt (or Peru) Current and its western 
continuation, the South Equatorial Current, is recognizable in the cloud pattern by a distinct zone 
of clear skies stretching from the Peruvian coast westward across the Galapagos Archipelago and 
along the equator. This clear zone seemingly follows the historical course of the Humboldt 
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Figure 1— ATS-II photograph taken on April 10, 1967, 1855 GMT, over the eastern Pacific Ocean. South 
America in the lower right is outlined by heavy cloudiness. The Peruvian coast is under clear skies; the 
Humboldt Current is marked by the extended band of no or little cloudiness through the center of the 
picture. Central America, the Gulf Coast, and southern Baja California are visible in the upper center 
portions. 


Current, which is a pronounced cold water stream that departs from the South American coast at 
Punta Negra, Peru. The current is, on the average, 2 to 5°K colder than the surrounding sea 
regions. This temperature difference seems to be quite sufficient to suppress the formation of 
convective cloudiness which is, in this case, present over all the surrounding warm water areas. 
In some portions, the boundaries of the cold water are even delineated by pronounced cloud bands 
that suggest the existence of a superimposed meso-scale thermal circulation between the cold 
and warm water regions, similar to the land-sea-breeze phenomenon. 
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In the series of photographs provided since November 1967 by the ATS-m satellite, which 
was stationary over the equator near South America, the outlines of the Humboldt Current are a 
pronounced phenomenon in almost every picture. Two examples (Figures 2 and 3) were selected 
in order to demonstrate the persistence of the exhibited cloud patterns from day to day. In both 
figures, the cold water of the Humboldt Current obviously suppresses the extended convective 
cloudiness present elsewhere; again, distinct cloud bands representing the upward branch of the 
induced local circulation accentuate at least parts of the current boundary. 

Depending on the season and the synoptic weather situation, the cold water area is not always 
free of clouds. Middle and high-layer cloud systems are influenced less by the thermal differ- 
ences at the underlying surface than are low-level and convective cloud systems. However, even 
in many cloudy cases, structural changes in the cloud patterns are very obvious over the Humboldt 
Current and help depict the location of the boundaries as is illustrated by Figure 4. In this case, 
large parts of the Humboldt Current are covered by mostly stratified clouds; nevertheless, the 
western current boundary can be recognized easily by a distinct difference in cloud density and 
even by a narrow break along the current boundary. This break can be understood as a reflection 



Figure 2 — Section of ATS— 1 1 1 photograph (IDCS) taken on March 11, 1968, 1622 GMT,show- 
ing South America and the eastern South Pacific Ocean. The largest parts of the continent 
are cloud-covered; the coasts of Chile and Peru are visible under clear skies that extend in- 
land from the clear region of the Humboldt Current. 
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Figure 3— Section of an ATS-lll photograph (IDCS) taken on March 12, 1968, 1659 GMT, 
showing South America and the eastern South Pacific Ocean. The Humboldt Current is de- 
marked by the extended band of clear skies through the center of the picture. 


of the subsiding branch of the local atmospheric circulation across the sea surface temperature 
discontinuity. 

From the foregoing discussion, it can be concluded that present satellite photography can pro- 
vide information useful in detecting pronounced sea surface temperature discontinuities, essen- 
tially along contrasting ocean surface currents. The location must be inferred visually from char- 
acteristic modifications of existing cloud patterns as caused by differential sea surface temperature 
patterns and the concurrently induced local vertical circulations in the lower atmosphere. How- 
ever, since the cloud features merely reflect these induced lower tropospheric local circulations 
rather than the sea surface temperature itself, the locations of the ocean current boundary and the 
cloud boundary need not coincide because of the width of the horizontal air motion across the cur- 
rent boundary. Therefore, an allowance for a displacement of approximately 10 kilometers toward 
the warmer water seems to be adequate. Also, a larger-scale atmospheric flow pattern might be 
superimposed on the local circulation, thus resulting in an additional displacement of the cloud 
boundary in a direction dependent on the superimposed wind field. The average distance between 
a distinct cloud boundary parallel to the Gulf Stream north wall from the sea surface temperature 
discontinuity, observed near the center of a la~ge atmospheric high pressure system and reported 
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earlier by Warnecke, et al. (1968), was approx- 
imately 20 to 40 kilometers. Additionally, mov- 
ies from ATS-HE time-lapse photographs taken 
over full daytime periods indicate a slight, but 
pronounced, diurnal variation of the cold pat- 
terns over the Humboldt Current or along its 
boundary. These effects may complicate the 
locating of an ocean current boundary accurately 
by its influence on cloud patterns, but no de- 
tailed and quantitative investigation of this 
problem has been accomplished. 


THEORY OF SEA SURFACE TEMPERATURE 
MEASUREMENTS FROM SPACE 
IN THE INFRARED 

Interpretation of the Radiative 
Transfer Equation 

The thermal radiation leaving the earth’s 
atmosphere at satellite altitude and measured 
through a filter system can be described by the 
radiative transfer equation 


N 


i r 2 

-J •(\)r i (\)e.(\)B(\. T.)d\ 

j f X 2 f h ® dr(K, h) 

+ - J <b(\) — B[A., T(h)] dh d\,(D 



Figure 4— Section of an ATS— 1 1 1 photograph (MSSC cam- 
era) taken on May 31, 1968, 1614 GMT, showing South 
America and the South Pacific Ocean. The coastal 
regions of Chile and Peru are visible under clear skies. 
The southern parts of the Humboldt Current are charac- 
terized in this picture by clear skies, while clouds cover 
the northern parts of the ocean current. The underlying 
sea surface temperature distribution is, however, re- 
flected in the cloud structure over the cold water 
current. 


where 

N is the detected radiance, 

\ is the wavelength, 

$(X) is the filter function = 0 outside the interval \ 1? \ 2 ], 

r(\) is the transmissivity from a height h to the satellite at height h*, 
r s (\) is the transmissivity from the radiating surface to the satellite, 
e s (\) is the emissivity of the radiating surface, 

B(\, T) is the Planck function, 

T is the absolute temperature, and 

T s is the temperature of the radiating surface. 
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The first term on the right side of this equation describes the contribution of an opaque sur- 
face, both solid and liquid (such as the earth, ocean, thick clouds, etc.), to the recorded radiance; 
the second term describes the radiation originating from the gaseous compounds of the atmosphere 
and optically tenuous constituents such as thin clouds, and dust. 

In an ideal atmospheric "window" (transmissivity t = t b = 1) Equation 1 reduces to 

i r ^ 2 

n w = I (2) 

x i 

With <E(\) known, the measured radiance from a water surface N w depends only on e s (X.) and 
T s . The emissivity of water is well known and varies little (between 0.90 and 0.99, depending on 
wavelength and view angle) in the presently-used infrared windows (Bramson et al., 1964). For 
most meteorological purposes, such as detecting clouds, black body radiation (e = 1.0) can be as- 
sumed without introducing a large error. Thus, N w depends on only one unknown yariable— i.e., the 
temperature T s of the radiating surface. In the case of clear skies and over water surfaces, T s 
closely represents approximately the sea surface temperature. If absolute temperature values, 
accurate within a few degrees Kelvin are required, corrections for e < 1.0 become necessary, at 
least on the basis of an average emissivity. 

Definition of Sea Surface 

On the basis of Kirchhoff T s Law, the high emissivity of liquid water in the infrared also means 
that its absorptivity is higb. A layer of less than 1 mm of water is opaque for infrared radiation from 
8 to 14 microns; 98 percent of the absorption occurs even in the first 0.1-mm layer (Ewing and 
McAlister, 1960). This me,ans that the observed infrared radiation emerges from the very top of 
the air-water interface, the "microsurface" where transfer processes (evaporation, reflection) 
take place and also microsurface properties (protein monolayer, dust, and oil films) are signifi- 
cant. Conventionally, sea surface temperature is defined as a mean temperature of a well-mixed 
layer at least several decimeters deep in reasonable accordance with the nature of conventional 
sea surface temperature observations ("bucket" or "intake" thermometric measurements). How- 
ever, experiments have shown (Ewing and McAlister, 1960) that the microsurface is approximately 
0.6°K colder than the layer immediately underneath (15-cm depth). It has also been observed that 
the microsurface is remarkably stable; it regenerates within tens of seconds after being disturbed 
by intense stirring or breaking of waves. Thus, under clear-sky conditions, radiometric meas- 
urements of sea surface temperature generally should be approximately 0.6°K lower than con- 
ventional measurements. 

Atmospheric Attenuation 

In reality, no atmospheric "window" is absolutely clean. Therefore, exact considerations 
must account for absorption of infrared radiation at least by the predominant tri-atomic atmos- 
pheric gases— water vapor, ozone, and carbon dioxide; then the complete radiative transfer equa- 
tion (Equation 1) will apply. According to Equation 1, absorption (and reradiation) vary with the 
gas distribution dr/dh and the gas temperature T(h) along the optical path; also absorption will vary 
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greatly with wavelength, depending on the actual location and extension of the absorption bands 
within the chosen filter range. Also, depending on the spatial distributions of the individual gases, 
their contributions to the detected outgoing radiance at satellite altitudes will peak in different 
layers in the atmosphere. Water vapor is the most important absorber and the most difficult one 
to handle because of its large variability in space and time in the atmosphere. Nevertheless, 
proper corrections must be applied to the satellite measurements in order to derive correct sea 
surface temperatures. These corrections will be discussed later in conjunction with the individual 
experiments. 


Cloud Interference 

A unique interpretation of the radiance measurements from a satellite in terms of sea surface 
temperature involves the assumption of a homogeneous, nontransparent target that fills the ra- 
diometer field of view. Because of the generally observed, vertical temperature decrease in the 
troposphere, any cloud "seen" by the radiometer will reduce the detected radiance and result in 
lower target temperature. This is true for opaque clouds filling a part of the field of view as well 
as for partly transparent clouds such as thin cirrus. The probability of having the field of view 
undisturbed by clouds increases with a decrease in the size of the field of view. Reducing the field 
of view, however, is limited by the sensitivity of the sensor because the signal strength decreases 
with a decreasing field of view; thus the signal-to-noise ratio becomes more and more unfavorable. 
Also the data volume increases with a decreasing field of view if total area I coverage is main- 
tained, thus imposing greater data reduction requirements. In order to determine whether clouds 
contaminate the radiation measurements, simultaneous cloud photographs from the same satellite 
will be a valuable source of information for proper interpretation of the radiation data. 


Daytime Observations 

The foregoing considerations of radiative transfer hold for thermal radiation only. These con- 
ditions are fully realistic, however, only on the night side of the earth. During daytime, they are 
still valid for infrared wavelengths beyond 5 microns. In the near infrared, 0.8 micron to 4 
microns, reflected sunlight is superimposed on the thermal radiation, and the term 


i 

R s (\) = - ®(\) r, (\)H(\)cos zd\ 


must be added on the right side of the radiative transfer equation (Equation 1). This term repre- 
sents the reflected solar radiation seen in the field of view and through the radiometer filter $(\) ; 
H(\) is the spectral solar irradiance; H(\) • cos z the spectral irradiance on the horizontal reflecting 
surface, where z is the zenith angle of the sun; and r s is the reflectivity of the viewed surface. 
Isotropic reflection is assumed. 
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In a perfect atmospheric window (r = 1,0), the radiative transfer equation then will apply in 
the form 

1 f^ 2 1 f^ 2 

N w = e s (^)B(X, T s ) d\ + - J $(\)r s (\)H(M coszdX. 

The radiometer, however, is calibrated in accordance with the thermal radiation of a black body 
of known temperature T: 



Thus, for a natural target, the simple assumption of T BB representing T s is no longer valid since 
with r s (\) a new variable is introduced. Also, taking into account the residual atmospheric absorp- 
tion within the window, H(\) will also become variable because of the variations in atmospheric 
transmission. Under these conditions, deriving sea surface temperature accurately from daytime 
radiometric measurements under clear sky conditions becomes very difficult. However, it will be 
shown later that certain characteristics of sea surface temperature patterns such as sharp temper- 
ature discontinuities at ocean current boundaries are recognizable in the near-infrared daytime 
observations. 

INFRARED RADIATION EXPERIMENTS 

Medium Resolution Radiometer Measurements 

Radiometers sensitive in the 8- to 12-micron M window” region were flown on Tiros II, III, IV, 
and VII. The large spectral range of approximately 4 microns involved an appreciable attenuation 
of the measurements predominantly by water vapor, but also by atmospheric ozone (9.6-micron 
band). Because the actual vertical distribution of these gases is known only at certain times and at 
relatively few places over the globe (i.e., in the vicinity of radiosonde stations which are particularly 
rare in the oceanic areas), it was necessary to develop correction models (Wark et al., 1962) on 
the basis of average conditions. These models are satisfactory for most meteorological purposes; 
however, for oceanographic applications, the uncertainties involved in the correction models are 
too large to yield reliable information on sea surface temperature. Additionally the 5-degree field 
of view of the medium resolution infrared radiometer (MRIR) resulted in a scan-spot at the earth's 
surface at least 60 kilometers in diameter, which is too large for oceanographic requirements. 

High Resolution Radiometer Measurements 


For the Nimbus satellites, a high resolution infrared radiometer (HRIR) with a 0.5-degree 
field of view responding in the 3.5- to 4.1 -micron window was developed, and it became possible 
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for the first time (Allison and Kennedy, 1967) to derive certain sea surface temperature features 
from the Nimbus I measurements. However, the short lifetime of the Nimbus I experiment se- 
verely limited the volume of useful data. 

The Nimbus II satellite, with an HRIR lifetime of 6 months, furnished a large volume of high 
quality data that could be explored more intensively for oceanographic applicability. A number of 
examples and preliminary results have already been published (Nordberg and Samuelson, 1965; 
Warnecke et al., 1968; Wilkerson, 1967; National Council on Marine Resources and Engineering 
Development, 1967; La Violette and Chabot, 1968; Rao, 1968; and Greaves et al., 1968). This re- 
port summarizes the results of a more systematic investigation. 


The Nimbus II HRIR Experiment 

The Nimbus II satellite was launched on May 15, 1966 into an earth- oriented, quasi-polar, near- 
circular, sun- synchronous orbit. Its perigee was 1095 km, the apogee 1179 km, and the orbital in- 
clination angle was 100.3 degrees. The ascending nodes occurred 32 minutes before noon, local 
mean solar time, and varied in position from day to day; graphs of the orbits are published in 
catalog form (Nimbus Project, 1966b). 

The high resolution infrared radiometer (HRIR) experiment lasted from May 20 through Novem- 
ber 15, 1966. The instrument, described in more detail in the Nimbus H User's Guide (Nimbus 
Project, 1966a), was a single -channel scanning radiometer sensitive in the 3.5- to 4.1-micron 
window region. Its sensor was a lead selenide (PbSe) photoconductive cell radiatively cooled to 
approximately 198°K. The radiometer was designed to measure radiance temperatures between 
210°K and 330°K with noise- equivalent-temperature differences of 1°K for a 250°K background. In- 
flight calibration was achieved by viewing alternately space and the spacecraft housing, the tem- 
perature of which was monitored during each scan. 

Ground calibration of the radiometer in terms of effective radiance 



was performed against a black body of variable, known temperature. The effective radiance meas- 
ured from space, 


N eff = f 4(\)N(\)d\ , 


thus is translatable into "equivalent black body temperature" T BB of a radiating source filling the 
field of view of the sensor. 
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Scanning was performed by a mirror rotating perpendicular to the orbital plane at a rate of 
44.7 rpm. Because of the field of view, satellite altitude, and satellite motion, continuous mapping 
was accomplished with contiguous scanning near the subsatellite track but with increasing overlap 
toward the horizons. The aperture and the height resulted in a rectangular scan spot approximately 
8.6 by 9.7 kilometers at the subsatellite point. 


Data Display 

A continuous reading of the radiometer output in analog form was transmitted from the space- 
craft. These data are available in three forms: 

1. as a function of satellite time in analog form, as recorded at the ground station, 

2. in the form of a playback on photographic film with a superimposed, approximate geographic 
grid, and 

3. as a digitized version of the radiometer data documented on Nimbus meteorological radia- 
tion tapes (NMRT) that can be used readily for computer processing, particularly for the 
production of two-dimensional data display in the form of standard map projections. 

The analog data contain the total volume of information. The tremendous data volume and the 
difficulty of analyzing analog data limit the practical application to selected studies where extreme 
detail is required. 

The photographic display of the analog data is the most convenient one, and a survey catalog 
of all available filmstrips has been published (Nimbus Project, 1967). For proper location of the 
data, a computer-processed, approximate geographic grid was superimposed on each filmstrip. 

The accuracy of the grid is nominally within 2 degrees of great circle arc at the subsatellite point. 
At the sides of the filmstrip near the horizons, severe distortions occur because of the scanning 
geometry. This effect and also the fact that an absolute relation between film density and temper- 
ature cannot be derived for these contrast-enhanced filmstrips further limit their applicability. 

Digitization of the analog data in conjunction with the concurrent orbital and attitude informa- 
tion leads to the mentioned NMR tapes, which represent the most useful data documentation for 
quantitative work. Mapping programs exist that print the data in standard map projections. The 
advantages of this form of data display are the presentation of absolute values, the elimination of 
distortion, and the possibility of automatically composing measurements from consecutive orbits 
into quasi -synoptic areal maps. However, because of the scanning geometry, either a loss of detail 
will result from averaging in the center portions of each swath or data gaps will occur at some 
distance from the subsatellite point where data spacing may be larger than the grid interval. 

These disadvantages, on the other hand, can be eliminated at least partly by the selection of proper 
map scales. Experience has shown that for oceanographic purposes a map scale of 1:1 million ap- 
pears to be most favorable. Automatic contouring or the transformation of the numerical values 
into easily distinguished color schemes are well proven additional aids for easier recognition of 
sea surface temperature patterns, while conserving the quantitative information. 
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Theoretical Data Accuracy and Interpretation 

The 3- to 4-micron window is located in the near-infrared. In this spectral region, solar ra- 
diation reflected from the surface of the earth and clouds is superimposed on the thermal radiation 
seen by a downward-looking instrument under daylight conditions. Therefore, the physical rela- 
tionship between the measured radiance and the temperature of the radiating surface is no longer 
valid. Consequently, in contrast to measurements in the 8- to 12-micron window, daytime Nimbus 
HRIR measurements are not applicable for the determination of sea surface temperature, although 
qualitative information such as the location of sharp sea surface temperature discontinuities can 
be derived under clear-sky conditions. 



Figure 5— The 4-micron window region. The upper solid curve represents the solar irradiance of a 6000°K 
black body at the top of the atmosphere, and the lower solid curve an infrared prismatic solar spectrum (Shaw 
et a!., 1951). The darker, shaded area represents the solar energy loss, mainly caused by molecular absorp- 
tion, in the earth's atmosphere. The horizontal lines indicate approximately the regions of absorption due to 
the major absorbing gases. The vertical lines indicate the center of absorption band due to major con- 
stituents (Goody, 1964). The dashed curve is the effective spectral response of the HRIR.(after Kunde, 1965) 
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The theoretical relationship between the measured nighttime equivalent black body tempera- 
ture and surface temperature measured by the Nimbus high resolution infrared radiometer as de- 
scribed for the Nimbus I instrument by Kunde (1965) is also valid for the almost identical Nimbus 
H HRIR. 


As indicated in Figure 5, residual molecular absorption in the 3- to 4-micron region is caused 
by water vapor (H 2 0), carbon dioxide (C0 2 ), methane (CH 4 ), and nitrous oxide (N 2 0). The absorp- 
tion by the minor constituents CH 4 and N 2 0 is negligible compared to the effect of C0 2 and water 
vapor, the spectral distribution of which is quantitatively shown in Figure 6. 

The results of computations of water vapor and carbon dioxide absorption for a number of 
model atmospheres and based on radiative transfer theory are reproduced in Figure 7. Absorp- 
tion by CH 4 and N 2 0 has been neglected and the emission and scattering by aerosols have not been 
taken into account. However, the results can be considered very close to natural conditions. 


From Figure 7 it follows that for zenith angles of the measured beam of less than 50 degrees, 
even in an extreme moist atmosphere (4.7-cm precipitable water content), the necessary correc- 



Figure 7— The fractional contribution to the total out- 
going effective radiance from each atmospheric layer. 
The assumed conditions are noted in the upper lefthand 
i corner (after Kunde, 1965). 


tion to the measured equivalent black body tem- 
perature does not exceed 4°K. On the average, 
a correction of +1 to 2°K will be sufficient. For 
higher accuracy of absolute data, an additional 
correction is necessary because of the deviation 
of the emissivity of water from 1.0 and the re- 
sulting reflection of atmospheric radiation. 


Assuming emissivity data from Bramson 
(1964), Table 1, and estimating the effect from 
Kunde ! s 1965 graphs, a depression of equivalent 



Figure 6— Theoretical surface temperature minus HRIR 
equivalent black body temperature difference as a func- 
tion of zenith angle for three model atmospheres. A 
grey body surface emissivity of unity and clear sky con- 
ditions have been assumed (after Kunde, 1965). 
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black body temperature from true surface tern- Table 1 

perature by 0.7°K results, on the average, for „ . . 

Emissivity of Liquid Water as a Function of Zenith 

zenith angles less than 60 degrees. The com- Angle and Wavelength in the 3- to 4- Micron Window 
bined atmospheric and emissivity effect, which (after Bramson, 1964). 

amounts to 2.7°K on the average ranges from 
+1°K (vertical view through a dry, high-latitude 
winter atmosphere) to 6.4°K (viewing at a 60- 
degree zenith angle through a very moist trop- 
ical atmosphere). 

In summarizing, the theoretically derived 
total correction adds up to +3.8 ± 2.0°K, aver- 
age. The single components are listed in Table 
2. This is in good agreement with earlier ex- 
perimental results from high-flying aircraft, as 
reported by Goldberg et al. (1964). 

Instrument Noise 

From the original analog records of two 
single scans of the spacecraft housing chosen at 
random during the period of the experiment, an 
average systems noise was determined to be 
±4°K. Contained in this noise are two com- 
ponents of particular interest— a permanent os- 
cillatory component and a sporadic interference 
of reflected sunlight. 

Although heavily reduced in amplitude when 
compared with Nimbus I HRIR measurements 
(Fujita and Bandeen, 1965), a residual oscilla- 
tory noise occurred in the Nimbus II HRIR 
measurements, as described by McMillin (1969). 

The noise frequency was found to be 200 Hz; the amplitude was about 2 C K. Because of strong in- 
terference of this noise with the high-resolution temperature patterns required for oceanographic 
or geological applications of HRIR data, a mathematical filter was produced by McMillin to 
eliminate this noise component before map production from Nimbus meteorological radiation tapes. 


Table 2 

Theoretically Derived Average Corrections to Nimbus 
II HRIR Measurements of Equivalent Black Body Tem- 
peratures for Derivation of Sea Surface Temperature 
Under Clear- Sky Conditions. 


Cause of Correction 

Temperature 

(°K) 

Microsurface effect 

+0.6 

Minor constituent absorption 

+0.5 

Water vapor and C0 2 absorption 

+2.0 

Emissivity effect 

+0.7 


This system noise is not reduced significantly by the routine averaging involved in the pro- 
duction of a large-scale grid print map; for example, in a 1:1 million map, only one to five in- 
dividual scan spots are averaged per grid square. Thus the numerical filter must be applied. The 
effect of the filter is shown in Figure 8. However, if the number of points in the average ap- 
proaches 15 or 20, as for example in a 1:3 million map, the spatial average itself becomes an 
adequate filter for most uses. 

Besides the permanently present oscillatory noise, a second, sporadic noise component ex- 
isted whenever the sun was close to the horizon. Direct sunlight was kept away from the sensor 
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Figure 8 — HR I R temperature patterns (in °K) derived from Nimbus II measurements on August 17, 1966 over the South 
Atlantic Ocean. The left one of the hand-colored figures shows raw data, the right one filtered data where the 
oscillatory noise has been removed. Notice the 'cleaner 1 pattern in the filtered version which also exhibits a num- 
ber of more organized small-scale patterns than the unfiltered version. 


by sunshields, but, as the measurements show, a small amount of sunlight scattered at parts of the 
spacecraft must have entered the radiometer near satellite sunrise and sunset. This scattered sun- 
light increased noise as well as the radiometer output, i.e., in measured equivalent black body tem- 
perature. The temperature increased approximately 6°K, average, as can be seen from Figure 9, 
in which the oscillatory noise was not eliminated. The scattered sunlight effect is visible in the 
upper portion of Figure 9, where the overall temperature increase of 6°K is clearly shown. In 
photofacsimile displays (Figure 10), the effect is noticeable because of a number of parallel, darker 
stripes across the affected portion of the picture. 

Thus, in the areas of satellite day-night transition, which actually represent two very small 
sections of each orbit, an additional correction must be applied in order to derive accurate quanti- 
tative results. 
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Figure 9— Digital color print of computer grid print map of Nimbus II HR!R measure- 
ments from orbit 1197 on August 13, 1966 over the Lake Michigan area. The zone of 
reflected sunlight interference during the satellite day-night transition is accentuated 
by the color selection. The predominantly red (land) and green (Lake Michigan) areas 
in the lower part of the picture are without interference. The predominantly blue 
(land) and Yellow-green (Lake Michigan and Lake Superior) areas to the north repre- 
sent the affected region. 
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Data Location 

In general, the grids that are automatically 
superimposed on the filmstrips are claimed to 
be accurate within 2 degrees of great circle arc 
at the subsatellite point. However, the grid lo- 
cation of the Nimbus II data was usually better, 
with errors exceeding 1 degree being rare. 

However, even an error of 1 degree at the sub- 
satellite point causes a large error in the dis- 
torted region at the side of the filmstrips. 

Contrary to the filmstrip grid, which used 
predicted orbital parameters to provide an ap- 
proximate, real-time rectification of the analog 

data, actual orbital parameters were used to Figure i 0 _Photofacsi m ile print of Nimbus II HRIR data 

produce the Nimbus meteorological radiation taken on August 13, 1966 on orbit 1 197 showing, as in- 

tapes (NMRT). Hence, the data location in the dicated, the area of reflected sunlight interference dur- 

ing satellite day-night transition, recognizable by a 
grid print maps produced from the NMRT is broad band of darker stripes across the filmstrip, 
claimed to be accurate within 0.5 degree of 

great circle arc. Although this was generally found to be true, a small uncertainty remains be- 
cause no corrections were applied for pitch and yaw; however, these were proven to be small. 
Roll errors were eliminated automatically by the nature of the data location procedure within each 
individual scan. During the digitization process, the center of the swath was determined from the 
occurrence of the horizon on either side; the center was then located at the subsatellite point 
(nadir). Some small residual uncertainty can still occur, however, because of uncertainties in the 
definition of the actual horizon. This results in a location error of ±0.5 grid intervals in a 1:1 
million map (1 grid interval corresponds to 1/8 degree longitude, i.e., 10 kilometers at 45 degrees 
latitude), as illustrated by Figure 11. In this grid print map showing equivalent black body tem- 
peratures around Lake Michigan close to local midnight on July 31, 1966, the contouring was 
chosen to accentuate the isotherm nearest the lakeshore. Because the eastern and western shore 
lines of Lake Michigan are fairly smooth and parallel, the small location error of the individual 
scan lines becomes recognizable. The data population in the map is 2 to 3 measurements, average, 
per grid point. Thus the general grid error of the entire map appears to be less than 0.2 degree 
of a great circle arc. 

SEA SURFACE TEMPERATURE PATTERNS DERIVED FROM NIMBUS II OBSERVATIONS 

A systematic search for sea surface temperature features in the Nimbus H HRIR data was 
performed for the most spectacular surface temperature discontinuity, the north wall of the Gulf 
Stream off the United States east coast. Nimbus II HRIR photofacsimile data for the 174 days of 
HRIR operation revealed at least parts of the Gulf Stream boundary in 48 nighttime and 9 daylight 
cases. It should be mentioned that the daytime number is not representative because daylight data 
were obtained only for selected days. Some of the 57 cases will be presented and discussed in 
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Figure 1 1 -Section of the computer grid print map of Nimbus I! HRIR nighttime measurements 
taken on July 31, 1966 (orbit 1024) over the Lake Michigan area. There are eight grid in- 
tervals per degree longitude. The superimposed grid represents the geographic grid produced 
by the computer for data location. Its accuracy is better than 0.2 degree of great circle 
arc. The outlines of the lakeshore are reflected by a sharp gradient in the printed equiva- 
lent blackbody temperatures (°K) as indicated by the drawn isotherms. The slight oscilla- 
tion of these isolines is caused by small location errors of the individual scan lines. 
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Table 3 


Indications of at Least Parts of the Gulf Stream Boundary from Nimbus II HRIR Measurements. 


Datr (1966) 

Data Orbit 

Remarks 

Date (1966) 

Data Orbit 

Remarks 

May 30 

198 

Fair 

October 1 

1849 

Good 

June 2 

238 

Good 

October 7 

1929 

Good 

June 3 

251 

Good 

October 8 

1492 

Excellent 

June 4 

265 

Fair 

October 9 

1955 

Poor 

June 6 

291 

Fair 

October 12 

1995 

Excellent 

June 7 

305 

Fair 

October 13 

2009 

Fair 

June 15 

411 

Good 

October 14 

2022 

Good 

June 19 

463 

Fair 

October 15 

2035 

Fair 

June 22 

504 

Fair 

October 16 

2048 

Fair 

June 28 

584 

Fair 

October 22 

2128 

Fair 

July 1 

634 

Fair 

October 28 

2208 

Fair 

July 4 

664 

Excellent 

October 29 

2221 

Fair 

July 17 

837 

Poor 

October 29 

2222 

Fair 

July 22 

904 

Poor 

October 30 

2235 

Fair 

July 23 

917 

Good 

October 31 

2248 

Poor 

August 1 

1037 

Good 

November 5 

2315 

Excellent 

August 2 

1050 

Good 

November 8 

2355 

Poor 

August 12 

1183 

Good 

November 10 

2381 

Poor 

August 19 

1276 

Fair 

November 15 

2448 

Fair 

August 28 

1396 

Good 




August 30 

1423 

Fair 

May 18 

44 

Fair, daylight 

September 3 

1476 

Fair 

May 18 

45 

Good, daylight 

September 12 

1596 

Fair 

May 20 

71 

Good, daylight 

September 13 

1609 

Poor 

May 23 

111 

Good, daylight 

September 23 

1742 

Poor 

May 24 

124 

Good, daylight 

September 24 

1756 

Excellent 

May 25 

138 

Fair, daylight 

September 27 

1796 

Fair 

June 22 

510 

Good, daylight 

September 27 

1795 

Fair 

June 22 

511 

Good, daylight 

September 30 

1836 

Fair 

June 24 

537 

Excellent, daylight 


this report; a complete listing of available Gulf Stream cases is given in Table 3. The remarks 
refer to how well the Gulf Stream boundary can be distinguished in the HRIR data. 


Daytime Infrared Observations 

Nine daylight cases indicating sea surface temperature characteristics in the Gulf Stream 
area are contained in the Nimbus II HRIR data volume (Nimbus Project, 1967). One of these and 
two other cases of indications of the Brazil and Falkland Currents boundary were selected for il- 
lustration in this paper. 


The Brazil and Falkland Current Boundary 

Photofacsimile pictures of the HRER measurements and television photographs (AVCS) taken 
near local noon on May 25, 1966 by Nimbus II are reproduced in Figures 12 and 13. In the HRIR 
picture, the light-to- darker grey tones indicate low-to-higher radiances. The combined temperature 
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Figure 12— Photo facsimile and schematic analysis of Nimbus II HRIR daytime measurements taken on May 25, 1966, 

1546 GMT (orbit 137), over the South Atlantic Ocean. 
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Figure 13— Mosaic and schematic analysis of Nimbus II AVCS television photographs taken on May 25, 1966 

1546 GMT (orbit 137), over the South Atlantic Ocean. 





and albedo contrast between land and water is clearly pronounced in the images of outgoing radi- 
ance, revealing fine details such as the Rio de La Plata, the Parana River system, the mouth of the 
Uruguay River, and the large lagoons near the coast of Southern Brazil. An interesting detail that 
was visible in the original photograph is that large cities such as Montevideo and Buenos Aires 
were indicated as dark spots of higher radiance because of higher surface temperature and higher 
reflectivity. Over the ocean, a sharp discontinuity in radiance was found along the boundary be- 
tween the southward-flowing, warm Brazil Current east of 50° to 53 °W and the colder Falkland 
Current heading north closer to the South American coast (Figure 12). A mosaic of the concurrent 
cloud photographs is given in Figure 13. By comparing simultaneous television and infrared ob- 
servations, cloudfree areas can definitely be depicted. Large parts of the discontinuity in the HRIR 
measurements were located under clear skies and therefore can be interpreted as sea surface 
temperature discontinuity although quantitative temperature values cannot be derived for reasons 
explained previously. The HRIR data in Figure 12 show that clouds can appear darker or lighter 
than surface features, depending on their actual temperature. The large cloud system consisting 
of two major cloud masses that extend southward from 40°S, 50°W on both sides of 50°W has the 
eastern cloud mass darker than the western, indicating that the eastern clouds are at a lower level. 
However, in the television pictures, Figure 13, both of these cloud masses appear almost equally 
bright. It is also seen in the HRIR data that the eastern cloud mass is also darker (meaning 
warmer) than the cold Falkland Current farther to the north. 

A phenomenon of particular interest in this case is the distinct cloud band about 100 km to 200 
km off the Brazilian coast between 26°S and almost 33 °S. This pronounced line of clouds is very 
near the temperature discontinuity (current boundary) in the HRIR measurements and seems to in- 
dicate, in the way described during the discussion of cloud photographs taken over the Humboldt 
Current, the upward branch of the local vertical motion system induced by differential heating 
along the line of sea surface temperature contrast. It should be considered, however, that the up- 
ward motion most likely occurs several kilometers inside the area of warmer water rather than 
at the discontinuity itself. This example demonstrates the advantage of simultaneous multispectral 
observations in complementing each other and thus favoring a more valid and useful interpretation 
of each of the observations. 

The second selected case shows the Falkland and Brazil Current boundary 1 month later, near 
local noon on June 24, 1966. In this case, the current boundary is not delineated by a cloud line 
and appears, for longer distances, under clear skies which favors the interpretation of the day- 
time HRIR observations (Figure 14). A discontinuity in radiance (grey tone) can be followed along 
a nearly meridional line across the entire picture, as indicated in Figure 14. Comparing this line 
with the cloud photograph in Figure 15 shows that the section between 32° and 39°S of this discon- 
tinuity is in clear air and therefore must represent a sea surface temperature feature. Farther 
south, this conclusion is no longer certain because the discontinuity follows the boundary of an ex- 
tended cloud system too closely. It is likely, although not conclusive, that this apparently low- 
level, convective cloud system is caused by slightly decreased vertical stability over the warmer 
water surface in an air mass of nearly adiabatic lapse rate. Such a case is theoretically plausible, 
and a very pronounced real case has been reported by Warnecke et al., 1968. However, no 


21 




60 °W 50 °W 60 °W 50°W 

Figure 14— Pholofacsimil e and schematic analysis of Nimbus II HRIR daytime measurements taken 
on June 24, 1966, 1500 GMT (orbit 536), over the South Atlantic Ocean. 


substantiating evidence can be brought in this case; therefore the current boundary indicated in 
Figures 14 and 15 must be considered uncertain south of the bar marked near 39°S. 


The Gulf Stream 

The third example of daytime HRIR data was taken over the Gulf Stream area also on June 24, 
1966. A photofacsimile of these measurements is shown in Figure 16. Again, the land and water 
contrast caused by the temperature and reflectivity contrast under clear -sky conditions is very 
obvious in the observed radiances. Geographic outlines of the United States east coast can be 
recognized easily (Figure 16). Over the ocean, the pronounced discontinuity extending northeast- 
ward from about 70 km off Cape Hatteras can be interpreted in most portions as the sea surface 
temperature discontinuity along the north wall of the Gulf Stream. Patches of cirrus, which are 
white in Figure 16 because of their low temperature and are easily identified by their structure 
and thin appearance in Figure 17, are located in only a few minor areas near or over the boundary 
to the extent that they mask its location. Toward the southwest, along the Carolina coast, the cur- 
rent boundary is not defined as sharply but can still be traced. The southeastern boundary of the 
Gulf Stream is naturally a much less defined feature in the HRIR picture, but at least indications 
can be seen of this boundary zone which seems to follow parallel to the meandering north wall. 
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Figure 15— Mosaic and schematic analysis of Nimbus II AVCS television photographs taken on June 24, 1966, 

1500 GMT (orbit 536), over the South Atlantic Ocean. 
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Figure 16— Photofacsimile and schematic analysis of Nimbus II HRIR daytime measurements taken on June 24, 1966, 

1528 GMT (orbit 538), over the North Atlantic Ocean. 
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Figure 17— Mosaic and schematic analysis of Nimbus || AVCS television photographs taken on June 2, 1966 

1528 GMT (orbit 538), over the North Atlantic Ocean. 

These three cases show that infrared measurements in the 3- to 4-micron window allow the 
observation of oceanographic features under daylight conditions from a satellite. In conjunction 
with simultaneous cloud photography, it is possible to outline the location of sharp sea surface 
temperature contrasts such as those existing along the north wall of the Gulf Stream and along the 
divide between the Brazil and Falkland Currents. 


Nighttime Infrared Observations 

Definition of the Gulf Stream Boundary 

The inshore boundary of the Gulf Stream was defined by Fuglister and Worthington (1951) as 
the line of no horizontal pressure gradient. Generally, the horizontal temperature discontinuity 
at 200-meter depth is considered as the actual current boundary location. The boundary surface, 
however, is inclined by approximately 1/3 to 1/2 degree from the horizontal plane. Thus, at the 
immediate sea surface, the boundary appears approximately 25 to 40 km north of the actual stream 
boundary (Stommel et. al., 1953, Stommel, 1965). This surface temperature discontinuity was ob- 
served by the remote sensing techniques and will be discussed here. 


The Gulf Stream North Wall 

Nighttime HRIR sea surface temperatures recorded over the northern edge of the Gulf Stream 
on June 2, 1966 are presented in Figure 18. This black-and-white filmstrip is based on the origi- 
nal data received at the ground station and therefore contains the oscillatory noise described 
earlier. The first-hand analysis shown in Figure 19 is based on a digitized version of the same 
data and therefore contains this noise component also. Because of the smoothing process involved 
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in this isothermal analysis, part of the noise 
seems to be eliminated; however, other parts 
might have been emphasized by this subjective 
process. In spite of this fact, the main charac- 
teristics of the sea surface temperature pattern 
stand out clearly and demonstrate the capability 
of the method. Although cloud interference by 
a cold front was observed over the main body 
of the Gulf Stream and sporadically in other 
parts of the figure, the north wall of the ocean 
current was well pronounced over a distance of 
approximately 700 kilometers. It is character- 
ized by a small zone of strong temperature - 
contrast of approximately 10°K, average, across 
a distance of 10 kilometers between 71° and 
76°W, as indicated by the approximately parallel 
heavy lines in Figure 19. This zone appears as 
a sharp discontinuity in grey tone in the photo- 
graph (Figure 18). The actual temperature in 
these infrared data was about 300°K over the 
warmer Gulf Stream waters, which appear 
darker in the picture, and about 290°K over the 
cooler slope waters to the north, appearing 
lighter in the photograph. 

A reanalysis of this case after use of a filter program to eliminate the periodic noise is re- 
produced in Figure 20. Although details changed, the overall picture remained the same, and the 
Gulf Stream boundary stood out even more clearly than before. Figure 21 shows the same data 
in an improved form of data display— the computer grid-print map in a Mercator projection of 
eight grid points per degree longitude with the third dimension (temperature) transformed into a 
selected color scheme. Each color square corresponds to the temperature at one grid point; the 
temperatures are translated into distinct colors in 2°K intervals. The Gulf Stream boundary, ac- 
centuated by the specific color selection used, is visible along the northern edge of the red and 
purple areas (warmer than 296°K) where a sharp transition to green (286° to 292°K) can be ob- 
served. With regard to the absolute values of the measurements presented in this case, it is im- 
portant to notice that from the northern margin southward to approximately Cape Hatter as inter- 
ference by sunlight reflected from spacecraft parts, as described earlier, occurred; it is indicated 
by the pattern of broad, darker streaks across the upper part of Figure 18 and by the shift of red 
to purple and light green to dark green in the upper part of Figure 21. 

The HRIR measurements on August 28, 1966 were selected as a second case of study. Al- 
though in this season the temperature contrast between the Gulf Stream and the slope water is 
significantly reduced, the inshore edge of the Gulf Stream is still pronounced in the photographic 
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Figure 18— Photofacsimile of Nimbus I! HRIR nighttime 
measurements on June 2, 1966, 0507 GMT (orbit 238), 
over the North Atlantic Ocean near the United States 
coast. 
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Figure 19— First analysis of Nimbus II HRIR nighttime measurements over the Gulf Stream taken on June 2, 1966 
(Allison et al., 1967), based on the unfiltered data. Isotherms are drawn in 10°K intervals; additionally, the 295°K 
isoline is shown. 
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Figure 20— Reanalysis of Nimbus II HRIR nighttime measurements over the Gulf Stream taken on June 2, 1966 (orbit 
238). Isotherms are drawn In 2°K intervals. The shaded area indicates those parts of the Gulf Stream seen by the 
radiometer without cloud contamination. 
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Figure 21— Digital color print of computer grid print map of Nimbus II HRIR measurements from orbit 238 on June 2 , 
1966 over the Gulf Stream (oscillatory noise eliminated by numerical filter). Blue and grey areas indicate clouds 
and/or land. The North Wall of the Gulf Stream stands out under clear skies; its southern boundary is obscured by 
clouds. Note the higher temperatures shown in Chesapeake Bay and the Delaware River if compared with the near- 
shore ocean temperatures. 


28 


I 


79 °W 69 °W 



79 °W 69 °W 

Figure 22— Photofacsimll e of Nimbus i! HRIR nighttime measurements taken on August 28, 1966, 

0440 GMT (orbit 1396), over the Gulf Stream area. The denumeration of the geographic grid 
takes into account an obvious error in the automatic picture grid. 

image (Figure 22) and in the colored version of the digital computer printout (Figure 23). The 
striking phenomenon of this case is the strong meander of the Gulf Stream between 67° and 73 °W, 
which is in good agreement with the documented Gulf Stream positions for this month (U. S. Naval 
Oceanographic Office, 1966). The Gulf Stream temperatures are shown between 296° and 300°K on 
the average, with a few spots warmer than 300°K. The slope water temperature ranges mainly 
from 292° to 296°K. This example demonstrates that even in the season with minimum tempera- 
ture contrast the HRIR method is capable of detecting the northern Gulf Stream boundary. 
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Figure 23 — Digital color print of computer grid print map of Nimbus II HRIR measurements from orbit 1396 on August 
28, 1966, 0440 GMT, over the Gulf Stream. The northern boundary of the Gulf Stream is indicated by the transi- 
tion from darker red (more than 296°K) to orange and yellow (292° to 296°K). Green, blue, and grey represent 
clouds and land. 


Comparison with Aircraft and Shipboard Measurements 

In both of the preceding cases, simultaneous radiometric observations of the sea surface tem- 
perature structure from low-flying aircraft were available for comparison. The airplane, oper- 
ated by the U. S. Naval Oceanographic Office, routinely maps the northern edge of the Gulf Stream 
between 65°W and 75°W. It is equipped with a Barnes infrared radiometer that is sensitive from 
7.3 to 13.5 microns, has a 2.2-degree square, fixed field of view, and in flight is accurate to ±0.4°K 
95 percent of the time. The airplane track is determined by Loran A. The navigational accuracy 
of Loran A over the Gulf Stream area is better than 4 km between 63 °W and 75 °N, according to the 
U. S. Naval Oceanographic Office (1967). 

The airplane usually performs two tasks. It tracks the inshore surface boundary of the Gulf 
Stream and records continuously the sea surface temperature profile across the detected boundary 
and over the other ocean areas to and from the scene. During the tracking portion of the mission, 
the aircraft is steered back and forth across the zone of strong temperature gradient with the aid 
of an infrared radiometer indicator in the cockpit. 
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The results of the research flight on June 2, 1966 are shown in Figure 24. The flight, con- 
ducted at night, commenced just before and extended slightly beyond the actual time of the Nimbus 
II overpass. The heavy dots are positions where the strong temperature gradient was located, and 
the solid lines with tick marks and temperatures are the tracks along which temperature profiles 
were taken. The temperatures shown are average values between two tick marks. These airplane 
data were superimposed on the almost simultaneously taken Nimbus II HRIR data analyzed from a 
high resolution computer printout. It can be derived from this figure that the zone of sharp tem- 
perature gradient in the HRIR data coincides surprisingly well with the locations detected by the 
aircraft, particularly if one considers the uncertainties involved in the location of both of the data 
sets. The average distance between the 27 marked Gulf Stream boundary locations from aircraft 
detection (black dots) and the center of the zone of strongest horizontal gradient of sea surface 
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Figure 24 — Comparison of simultaneous Nimbus II HRIR (orbit 238, 0500 GMT) and aircraft radiometric nighttime 
measurements (0300 to 0630 GMT) of sea surface temperature near the northern Gulf Stream boundary on June 2, 
1966. Heavy dots represent the location of the Gulf Stream boundary as detected from the aircraft. Numbers along 
the tracks of specially flown horizontal temperature profiles represent average sea surface temperatures (in °K) be- 
tween tick marks. 


31 














temperature detected from the satellite is 12 kilometers. For the temperature comparison, some 
data near the strongest gradient zone had to be excluded in order to diminish mere location errors. 
The average temperature difference between the aircraft and the satellite measurements for the 
remaining 94 pairs of simultaneous data was determined to be 0.45°K; that is, the satellite data 
were, on the average, about 0.5 degree too warm. From the theoretical considerations previously 
discussed, the satellite data should give values that are too low by approximately 4°K (see Table 2). 
Because of interference caused by sunlight scattered from parts of the spacecraft, an error of +6°K 
must be assumed in this part of the orbit. Thus, both systematic errors appear largely to compen- 
sate each other, and sea surface temperature should be indicated 2°K too warm by the HRIR meas- 
urements. The observed difference of +0.45°K can well be considered reasonably small and within 
the limits of accuracy of this error analysis. 

In the southern part of the Gulf Stream region of Figure 20, ship data were also available for 
comparison. Research vessel EASTWARD of the Research Triangle Institute, North Carolina, was 
cruising across the Gulf Stream between May 23 and 28, 1966. Its temperature measurements are 
plotted in Figure 25. With regards to the sea surface temperature distribution, these measurements 
are fairly representative for the day of the Nimbus II overpass, June 2, 1966. Unfortunately, most 
of the ship data were located in the area where the satellite data were disturbed by clouds. Thus, 
a direct comparison was possible in only a few locations. Nevertheless, the ship measurements 
verified that the Gulf Stream temperatures were 297° to 299°K, thus confirming the airplane and 
satellite measurements. The inner edge of the Gulf Stream was found near 34.6°N, 76°W by the 
EASTWARD, and the Nimbus n data, taken 10 days later, show the strongest horizontal temperature 
gradient in almost the same position. 
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Figure 25— Comparison of Nimbus II HRIR nighttime measurements on June 2, 1966 (orbit 238) with temperature 
measurements from shipboard (R/V EASTWARD) between May 23 and 28, 1966. 




Figure 26— Comparison of Nimbus II HRIR nighttime measurements on August 28, 1966, 0440 GMT (orbit 1396), with radiometric aircraft measurements 
on August 27, 1966, 1200 - 2100 GMT. Numbers along airplane tracks are temperatures (in °K) averaged between tick marks. 


On August 28, 1966, data from another research flight over the Gulf Stream were available 
for direct comparison with Nimbus II measurements. These flight data, unlike those of June 2, 
were taken during daylight hours approximately 6 to 18 hours after the Nimbus II overpass. In 
this case, no interference by scattered sunlight was observed in the satellite data. The air- 
plane track and the results of the temperature measurements taken along these tracks are re- 
produced in Figure 26. A numberical evaluation for 196 pairs of simultaneous measurements with- 
out apparent cloud contamination of the satellite data revealed that the Nimbus II HRIR measure- 
ments were 3.9°K colder than the concurrent airplane data. This temperature difference is within 
0.1°K of the theoretically expected value. However, the theoretical considerations were based on 
standard atmosphere conditions, and it is reasonable to assume that in different seasons and 
locations, larger differences between theoretical and observed data may occur unless more re- 
alistic atmospheres are used for the theoretical computations. Ideally, corrections were to be 
computed from simultaneously taken spectrometric measurements of the actual vertical water 
vapor profile as suggested by the Nimbus III IRIS (Infrared Interferometer Spectrometer 
experiment). 


The Brazil and Falkland Current Boundary 

In addition to the daylight cases presented earlier, two cases of nighttime observations— only 
24 hours apart— of the boundary between the warm Brazil and the cooler Falkland Current were 
selected from the Nimbus II HRIR data for presentation. On August 16, 1966 behind nor th eastward - 



Figure 27— Phofofacsimile of Nimbus II HRIR nighttime 
measurements on August 16, 1966, 0245 GMT (orbit 
1235), over the Atlantic coast of South America. 


moving cold front, large parts of this boundary 
were under clear skies. Besides other very 
prominent features such as the South American 
coast, the La Plata River mouth, and details 
such as the Parana and Uruguay Rivers, the 
Brazil Current boundary is clearly visible from 
3 2 °S to 38 °S in Figure 27, Farther south, scat- 
tered cloud decks obscured the current bound- 
ary. Temperatures of 289°K and 281°K can be 
derived from Figure 28 as average surface tem- 
peratures of the Brazil current and the Falkland 
Current, respectively. After addition of the cor- 
rections, discussed earlier, these temperatures 
are in good agreement with climatological aver- 
ages for August in this sea region (Hydrographic 
Office, 1954). 

On the following day, August 17, 1966, this 
ocean current boundary was seen by the infrared 
radiometer between 30°S and 44 °S (Figure 29), 
i.e., over a much greater area. The northern 
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Figure 28 — Digital color print of computer grid print map of Nimbus II HRIR measurements oh August 16, 1966, 0245 
GMT (orbit 1235), over the Atlantic coast of South America. Oscillatory noise was removed from the data. Colors 
are separated in steps of 2°K. 


parts, which are seen on both consecutive days, are a good example of the repeatability of the ob- 
servation method. 

The equivalent black body temperature over the Falkland Current in the lower part of the fig- 
ure is around 280°K; a steady, slight increase can be noticed toward the north in Figure 30, where 
282°K can be derived, on the average, for the northern part of the Falkland Current. The Brazil 
Current, on the other hand, is characterized mainly by 288° to 290°K. An average temperature con- 
trast of 5.4 °K over a mean distance of 24 kilometers can be derived along the inshore edge of the 
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Brazil Current between 32° and 40°S for this day. 
At certain locations, this difference is as high as 
10°K, according to Figure 30. Another example 
for the clear definition of the inshore edge of 
the Brazil Current is given in Figure 3 1 by the 
reproduction of an original analog record of a 
horizon-to-horizon scan across the current 
boundary on August 17, 1966. Instrumental noise 
is a very striking feature in this analog record, 
but the temperature contrasts visible in the photo- 
facsimile display (bottom of Figure 30) can be 
identified easily. The numerical filter (developed 
by McMillin, 1969 to remove the predominant 
200- Hz oscillatory noise component of the dig- 
itized data) was used on a section of the same 
scan and is presented in Figure 32. Here the 
current boundary clearly stands out as a tem- 
perature step of almost 10°K. 


50 °W 


Figure 29— Photofacsimile of Nimbus I! HR1R nighttime 
measurements on August 17, 1966, 0210 GMT (orbit 
1248), over the South Atlantic Ocean. 


These sporadic observations indicate that 
the Brazil and Falkland Current boundary, which 
is part of the Subtropical Convergence in the 
South Atlantic Ocean, is as sharply defined as 
the inshore edge of the Gulf Stream and represents one of the outstanding features in the sea sur- 
face temperature distribution. 


The Oyashio-Kuroshio Region 

The Oyashio-Kuroshio boundary was seen in the Nimbus II HRIR data in a number of cases. 
Measurements taken on October 6 and 7, 1966 over this region were selected for this presentation 
because of the characteristic thermal structure exhibited in the area of interference between the 
two ocean currents. Large parts of northern Japan and the surrounding seas had clear skies on 
October 7, 1966. Contrary to the examples presented earlier, a number of interesting features, 
not only in connection with the main current boundaries but also within the main water boundaries 
themselves, can be observed in the temperature pattern exhibited by the photographic display of the 
infrared measurements (Figure 33). Most prominent is the temperature contrast between the cold 
land and warmer ocean waters near midnight. Within the ocean areas, a striking temperature con- 
trast can be seen between the sea regions on both sides of the Kuril Chain. According to the season, 
the Okhotsk Sea is actually still warmer than the Oyashio flowing southwestwardly east of the Kuril 
Islands. Very obviously the chain of islands prevents a strong mixing of water bodies of different 
temperatures. Farther south, the transition from the cold Oyashio waters toward the warmer 
Kuroshio seems to take place in two distinct steps. A sharp line of temperature contrast can be 
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Figure 30— -Digital color print of computer grid print map of Nimbus II HRIR measurements on August 17, 1966, 0210 
GMT (orbit 1248), over the South Atlantic Ocean. Colors are separated in 2°K steps. Oscillatory noise has been 
removed from the data. 


seen near 41°N; a second one, which seems to represent the northern edge of the Kuroshio itself, 
is very pronounced somewhat south of 40°N. From Figure 34 it can be seen that the equivalent 
black body temperatures over the Kuroshio were near 298°K, in the transition zone around 292°K, 
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Figure 31— Example of an analog record of an individual scan line of Nimbus II HRIR nighttime 
(orbit 1248) data taken on August 17, 1966 at 02 h 1 l m 42 s GMT across Argentina and the South 
Atlantic Ocean. The photofacsimile underneath indicates the exact location of the scan line. 



Figure 32— Section of the Nimbus II HRIR analog record 
reproduced in Figure 31 after digitization and applica- 
tion of a numerical filter to remove oscil latory noise. 


and about 286 °K in the southern parts of the 
Oyashio. Because of the possibility of some 
cloud contamination by perhaps thin cirrus or 
small convection clouds not resolved by the 
HRIR instrument, the observations on the pre- 
ceding day were inspected for the reproduci- 
bility of the observed patterns. Figures 3 5 and 
36 show these observations. Although a stronger 
cloud interference existed on October 6, 1966, 


38 


several of the characteristic patterns seen on 
the next day could be identified. Figure 37 
schematically indicates the most prominent 
structural features repeated in both sets of 
data. From this, it is obvious that the two 
boundaries observed are really phenomena in 
the sea surface temperature distribution. No 
conclusions on the history or dynamics of this 
dual boundary between the Oyashio and the 
Kuroshio can be drawn from these observations 
yet, but infrared data taken 1 month later over 
the same region (Figure 38) indicate the per- 
sistence of a very complex thermal structure 
of the Oyashio-Kurashio interface. 

The Agulhas Current 


130° E 140°E 



130°E 


140°E 


Figure 33— Photofacsimil e of Nimbus II HR I R nighttime 
measurements on October 6, 1966, 1500 GMT (orbit 
1921), near Japan. 


The Agulhas Current is a warm-water cur- 
rent flowing southwestward along the east coast 
of South Africa and, according to textbooks 
(Defant, 1961; Pickard, 1964), splitting near 
33 °S into a weaker branch flowing around South 

Africa into the South Atlantic Ocean and a stronger branch turning south and merging into the 
southern Subtropical Convergence near 40°S. Sharp temperature discontinuities within this ocean 
current were reported by radiometric measurements from aircraft off Durban, South Africa 
(Harris and Stavropoulos, 1967), but no specifically sharp current boundaries such as those as- 
sociated with other currents are shown in standard sea surface current maps (Defant, 1961). 


The Nimbus II HRIR data contain a good example, with sufficiently reduced cloudiness, that 
permits the observation of this particular sea region. The infrared photograph of September 19, 
1966, reproduced in Figure 39, shows a distinct sea surface temperature discontinuity from about 
33.5°S, 27°E to approximately 38°S, 20°E. This line seems to represent the meandering western 
edge of the Agulhas Current. Over the warmer water, equivalent black body temperatures of 292° 
to 294 °K were recorded, corresponding to sea surface temperatures around 298°K (after correc- 
tion); the corresponding values over the colder water were 285°K (T BB ) and 291°K (corrected 
temperature). However, as the pictorial display indicates cloud bands were still contaminating 
the area to some extent; thus these temperatures may be biased toward lower values. These 
measurements agree satisfactorily with climatic data. 


This example proves that structural features of the Agulhas Current can be detected and 
studied by observations from satellite altitudes. The sharpness of the temperature contrast found 
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Figure 34— Digital color print of computer grid print map of Nimbus II HRIR nighttime measurements on October 6, 
1966, 1500 GMT (orbit 1921), near Japan. Oscillatory noise has been removed from these data. The color steps 
are equivalent to 2°K temperature intervals. 
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in this case suggests the existence of a more 
strongly defined western boundary of the 
Agulhas current than was expected before. The 
temperature contrast, indicated by the high 
resolution infrared measurements of Nimbus, 
at least in the present case, seems to be as 
significant as the temperature discontinuities 
found along the Gulf Stream and Brazil Current 
boundaries. 

CONCLUSIONS 

Results from various meteorological sat- 
ellites suggest many possibilities for using 
orbiting platforms for oceanographic purposes. 
This paper presented a survey on the presently 
available observation techniques and data in 
sensing sea surface parameters from satellite 
altitude. 
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Figure 35— Photofacsimll e of Nimbus II HRIR nighttime 
Inference of sea surface temperature pat- measurements on October 7, 1966, 1430 GMT (orbit 

terns from meteorological satellite cloud pho- near Japan# 

tography seems to be a promising field for 

further investigation. The examples given show an obvious correlation between certain sea sur- 
face temperature discontinuities and the induced cloud patterns; however, detailed and more 
quantitative studies in conjunction with ship or airplane observations for ground confirmation are 
necessary in order to establish a reliable analysis method under cloudy conditions. 


In the case of clear skies, infrared measurements of high spatial resolution in one of the at- 
mospheric windows present a second method of deriving sea surface structural features from 
space. In the 3- to 4-micron window, the Nimbus II HRIR experiment acquired a sizeable quantity 
of data showing sea surface temperature features, sensed under daytime and nighttime conditions. 
However, in this near-infrared window, reflected sunlight does not permit the derivation of quanti- 
tative sea surface temperature data during daytime. Nighttime observations furnish quantitative in- 
formation not only on the location of certain ocean current boundaries, but also on the sea surface 
temperature distribution. Therefore, it would be desirable to change to the 8- to 12-micron win- 
dow, which is not affected by solar radiation, but lies entirely within the spectral range of thermal 
radiation of the earth- atmosphere system. Also it would be desirable, under daytime conditions, 
to have simultaneous, supplemental information on the cloud distribution from a second radiometer 
channel in the visible spectrum, thus permitting a more accurate interpretation of the infrared 
measurements in terms of sea surface temperature. During nighttime, no supplemental cloud ob- 
servation is available yet, although observations in the visible airglow bands from the OGO-4 
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Figure 36— Digital color print of computer grid print map of Nimbus II HRIR nighttime measurements on October 7, 
1966, 1430 GMT (orbit 1934), near Japan. Oscillatory noise has been removed from these data. The color steps 
are equivalent to 2°K temperature intervals. 
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Figure 37— The most significant features in the Nimbus 
II HRIR nighttime data near Japan which were repeated 
on October 6 and 7, 1966, representing landmarks and 
certain temperature discontinuities in the Oyashio- 
Kuroshio mixing region. 


Figure 38 — Photofacsimil e of Nimbus II HRIR nighttime 
measurements on November 5, 1966, 1420 GMT (orbit 
2320), near Japan. 
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Figure 39 — Photofacsimile and schematic analysis of Nimbus II HRIR nighttime measurements on September 19, 1966, 

2130 GMT (orbit 1698), near South Africa. 









satellite (Warnecke et al., 1968) and the development of a HRIR system to be flown on the ATS-F 
satellite offer promising possibilities for future research. 

Before techniques in other spectral regions such as passive microwaves are explored fully 
and developed, quantitative information on sea surface temperature by infrared radiometry will 
be restricted to clear regions. However, fairly useful sea surface data coverage can be achieved 
if measurements of clear areas taken over a period of several (perhaps 5 to 7) days are composed. 
This approach should be possible because of the 12-hour frequency of measurements at a certain 
location by one near-earth satellite and the relatively fast motion of cloud systems compared with 
the normally observed high persistence of sea surface temperature patterns within periods of 
several days. 

The examples presented of remote detection of sea surface temperature and ocean current 
characteristics, taken from satellite observations which were originally designed for meteorolog- 
ical purposes only, demonstrate that existing techniques are suitable, after some minor modifi- 
cations, to establish a satellite -borne, global oceanographic observation system that will open a 
new dimension of oceanographic research and application. 
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